Spin squeezing, entanglement and coherence in two driven, dissipative, nonlinear 
cavities coupled with single and two-photon exchange 
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We investigate spin squeezing, quantum entanglement and second order coherence in two coupled, 
driven, dissipative, nonlinear cavities. We compare these quantum statistical properties for the 
cavities coupled with either single or two-photon exchange. Solving the quantum optical master 
equation of the system numerically in the steady state, we calculate the zero-time delay second- 
order correlation function for the coherence, genuine two-mode entanglement parameter, and an 
optimal spin squeezing inequality associated with particle entanglement. We identify regimes of 
distinct quantum statistical character depending on the relative strength of photon-exchange and 
nonlinearity. Moreover, we examine the effects of weak and strong drives on these quantum statistical 
regimes. 
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I. INTRODUCTION 

Multiphoton processes in quantum optical systems 
have been intensely studied [2, 0] due to their central 
role in modern applications such as quantum switch- 
ing [H, quantum communication and computation 
They are also used for fundamental explorations of phase 
transitions in coupled nonlinear cavity or superconduct- 
ing (SC)circuit quantum electrodynamics (QED) sys- 
tems d-Gl]. 

An earlier study of two nonlinear cavities coupled with 
single-photon exchange [2fj| | revealed a curious interplay 
between coherence and localization of the photons. It 
is concluded that photons are coherent and delocalized 
over the cavities when the tunneling exchange is stronger 
than the nonlinearity. In the opposite case of weaker 
tunneling, photons are localized in each cavity and anti- 
bunched [21| . 

In the present work, we address the question if such an 
interplay can go across to quantum correlations as well. 
We specifically ask how such mutual influences between 
localization, coherence and quantum entanglement [22[ 
change under two-photon exchange f23T - |25j as well as un- 
der strong drive conditions [26l. |27]. Realization of two- 
photon exch ang e coupling was proposed for circuit QED 
systems |28l - l33| . 

Profound relations between different definitions of spin 
squeezing and entanglement have already been explored 
in det ail pj -Hj. Using Lipkin-Meshkov-Glick (LMG) 
Model [39], interplay of quantum correlations in squeez- 
ing and entanglement in finite quantum systems was dis- 
cussed quite recently |35[. Entanglement dynamics and 
exact properties of LMG model in the thermodynamical 
limit have been carefully analyzed in Refs. [40f- 42] . Two 
photon exchange coupled cavities can also be described 
by the LMG model [24[ . Our present contribution estab- 



lishes further connections between coherence and mode 
entanglement to spin squeezing and pairwise entangle- 
ment in this model under drive and dissipation. 

We calculate and compare the zero time delay sec- 
ond order quantum coherence function [43| . an optimal 
spin squeezing inequality associated with pairwise en- 
tanglement [441 ] and the genuine mode entanglement pa- 
rameter [45|]. Mode entanglement occurs in the second 
quantization description of the system; and hence it is 
fundamentally different from the particle entanglement, 
which happens in the first quantization description [46r - 
HH . Some possible realizations and experiments to detect 
multiparticle entanglement via optimal spin squeezing in- 
equalities can be found in Refs. |49rl5ll |. 

This paper is organized as follows. In Sec. Ql] we de- 
scribe the single and two photon exchange coupled non- 
linear cavity models under consideration. In Sec. Mil we 
introduce the parameters that we calculate to charac- 
terize coherence, entanglement and spin squeezing and 
present their steady state results. We conclude in Sec.lIVI 



II. THE MODEL SYSTEM 

We consider a system of two identical nonlinear cavi- 
ties, labeled with i = 1,2, coupled either by single or two 
photon exchange interactions. Both cavities are driven 
by a coherent pump at rate F at the laser frequency wl ■ 
The corresponding model Hamiltonians in a frame rotat- 
ing at uil can be written as 

= J2(Ubmhihi + Fb\ + F%) 



i=l,2 



+ J(b\b2 + blh), 
= ^2(Ublbibibi + Fbl + F*bi) 



J{bTbl 



iTbi), 



(i) 



(2) 
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where U is the nonlinearity parameter and J is the pho- 
ton exchange coefficient. Here we denote the model with 



2 



single and two photon exchange coupling as and 
H^ 2 \ respectively. These models describe generic two- 
site Kerr-Hubbard type interactions that may be realized 
in settings other than coupled cavities, for example in ul- 
tracold atoms [52-56]. Annihilation (creation) operator 
for the cavity photons with frequency Ui is denoted by bi 
(hi). We take Ui — u> L . 

Let us re-express the model Hamiltonians using the 
pseudo-spin operators of the cavity fields, 



•3 x 



J ii 



J z 



-^(blh-blh), 

\{b\k-b\b 2 ). 



(3) 



They satisfy the SU(2) spin algebra [J a ,Jp] = e Q ^J 7 . 
Here a, /3,7 £ {x, y, z} and e Q ^ 7 is the Levi-Civita den- 
sity. We obtain 



= UJ 2 Z 

ij( 2 ) = UJt 



2JJ X + Fb\+F%, 
2J(Jl-P)+Fb\+F% 



(4) 
(5) 



While the tunneling term is a mere rotation operator in 
the case of single photon exchange, it is a generator of 
spin squeezing that redistributes the spin fluctuations by 
twisting them about the two axis 57| in the case of two 
photon exchange. The nonlinear term always acts as a 
generator of spin squeezing by twisting the fluctuations 
around a single (z) axis. Spin squeezing is associated with 
multi-particle entanglement. Another type of entangle- 
ment can be found between the cavity modes. This mode 
entanglement is enforced by the mode coupling character 
of the tunneling terms. Drive term brings coherence into 
the system. 

In the next section we define the measures of quan- 
tum correlations to calculate spin squeezing, entangle- 
ment and coherence properties of these models. All the 
measures are closely related to the spin noise and hence 
strong interplay between coherence, squeezing and types 
of entanglement is expected. 



cavities. T>[x]p^> = [2xpv>x' — x'xfr- J > — p( J 'x'x]/2 are 
the Liouvillian superoperators in the Lindblad form. We 
assume k% = n and coherent pump amplitude F is taken 
to be real. 

We solve the master equation for the steady state den- 
sity operator using the QuTiP: Quantum Toolbox in 
Python software [58j . The master equation is solved by 
taking Ni = N 2 = 4 for the Fock space dimensions of 
each cavity field. When we increase the Fock space di- 
mension up to 8 and remake our calculations, we find that 
the results remain unchanged. We consider two cases dis- 
tinguished by the action of weak F/n = 0.1 MHz and 
strong pump F/k = 1 MHz. We take k/2tt = 0.4 MHz. 
The parameters we use in the simulations arc within the 
ranges accessible in present circuit QED systems [59|, [6(| ■ 

Numerically computed steady state density operator 
is used for calculation of second order coherence, spin 
squeezing and mode entanglement parameters. Zero-time 
delay second order quantum coherence function is defined 
by Hi 



5 (2) (0) 



tT(ttttbbp) 
[tr(StSp)]2 ' 

1 , ((An) 2 ] 



(h)< 



(7) 
(8) 



Here h = b'b, Ah = h — (h); and we use p instead of 
pU) f or notational simplicity. We evaluate g^(0) for 
b = bi. Both cavity photons have identical coherence 
functions due to exchange symmetry of our models un- 
der b\ «-> b 2 . g^(0) < 1 implies antibunching and sub- 
Poissonian statistics of photons, while g^ 2 \0) > 1 im- 
plies bunching and super-Poissonian statistics. Poisso- 
nian statistics and coherent photons are recognized by 
(/ 2 )(0) = 1. Antibunching condition is a violation of 
the Cauchy-Schwartz inequality obeyed by classical light; 
and hence antibunching is a profound manifestation of 
quantum light. 

Spin squeezing is witnessed by the inequality [44j 



(Jl) 



(J?) 



N 



< (N-l)(AJ m ) 2 



(9) 



III. RESULTS AND DISCUSSIONS 

To investigate the quantum dynamics of our model sys- 
tems under dissipation, we assume that the coupling of 
the cavity photons and the reservoir photons is weak; and 
the correlation time of the reservoir photons is negligibly 
short. Under these so called Born and Markov conditions, 
the dynamics of the open systems can be determined by 
solving the master equations 



ft) = 



-i[HV\pV>] 



E 



KiT>[x]p 



(6) 



for both single (j = 1) and two photon (j = 2) exchange 
cases. Here, p^ is the density operator for the corre- 
sponding case; ft, are the photon loss rates out of the 



where k, I, m take all the possible permutations of x, y, z. 
Inequality given in Eq. §§§ is one of the four optimal spin 
squeezing inequalities introduced in Ref. 44| with N be- 
ing the total number of particles in the system. Violation 
of it implies spin squeezing and particle entanglement. 
We rewrite the optimal spin squeezing inequality given 
in Eq. (|9j) with k = x, I = z and to = y as 

C = (J 2 X ) + {Jl) ~j-(N- 1)(AJ,) 2 < 0. (10) 

The positive values of £ indicate spin squeezing and pair- 
wise entanglement in the first quantization description. 
As we have an open system, N = (N) = (hi) + {h 2 ), with 
N = hi + h 2 and hi = bl&j, is not conserved and varies 
in time. 
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FIG. 1: (Color Online) Dependence of (a) ff (2) (0), (b) C (c) A12 for F/k = 0.1 and (d) £ (2) (0), (e) ( (f) A12 for F/k = 1 with 
respect to dimensionless J/k and U/k in two-site KH system with single- photon exchange. £ is multiplied by 10 3 and A12 is 
multiplied by 10 5 for visibility. 



Pairwise entanglement correspondence of spin squeez- 
ing can be difficult to interpret for photons due to lack 
of a commonly accepted definition of photon wave func- 
tion in position representation, especially in free space. 
Here we simply assume that the localized cavity fields can 
be considered as photonic wave functions. As the pho- 
tons are identical particles, their pairwise entanglement 
may be ambiguously determined by the spin squeezing in- 
equality for some cases and moreover it may not be a use- 
ful resource for quantum metrology applications (48l . loll ] . 
In the spirit of Ref. [62J, however, this should not pre- 
vent us to calculate spin squeezing for the purpose of 
distinguishing pairwise entanglement from the mode en- 
tanglement and to explore their relation to coherence. 

We investigate the entanglement in the second quan- 
tized description in terms of the genuine two-mode en- 
tanglement parameter A12 defined by [4o| 



\i2 = \(b\b 2 )\ 2 -(nin 2 ), 



(11) 



where hi — b\bi with i = 1,2. Positivity of the parameter 
A12 indicates mode entanglement. 

An intuitive link between coherence, spin squeezing 
and mode entanglement can be formed by considering 
that coherence function is a relative measure of number 
fluctuations which carries information on deviation from 
Poisson statistics. One can establish an analogy between 
the number operator h and the z-component of the an- 
gular momentum operator J z using h\ = N/2 + J z . The 
coherence function is influenced by squeezing the spin 
noise around the z axis. Under stringent conditions spin 
squeezing implies particle entanglement. Particle like be- 
havior, enforced by the photon antibunching can be as- 



sociated with spin squeezing type pairwise nonlinear in- 
teractions. In that case coherence and spin squeezing 
interplay can be further extended to multiparticle entan- 
glement. Similarly A12 can be expressed in terms of spin 
fluctuations and thus would be influenced by the coher- 
ence and spin squeezing in the system. 

Above heuristic arguments motivates the existence of 
an interplay between coherence, spin squeezing, pairwise 
and mode entanglement. Let us now analyze specific 
cases numerically. We will first consider the single pho- 
ton exchange case, then proceed to two photon exchange 
case in the following subsections. 



A. Single-Photon Exchange 



In Fig. [T(a)| Fig. [1(b) 



and Fig. 1(c) we plot sec- 
ond order coherence, spin squeezing and mode entangle- 
ment parameters as functions of nonlinearity and pho- 
ton exchange coefficients in the case of weak drive with 
F/2ir = 0.04 MHz, respectively. Fig. 1(a) reproduces the 



result in Ref. [20j . Second order coherence varies over the 
range of 0.00183 < g {2} (0) < 1.011. Coherent derealiza- 
tion of cavity photons happens for J/k > 1 and J > U, 
which corresponds to the photon bunching and Poisso- 
nian statistics indicated by </ 2 '(0) ~ 1. Strong photon 
antibunching regime with # (2) (0) ~ lies in the region of 
U/k > 1 and J/k < 1. 



Fig. 1(b) shows that spin squeezing always present in 
the system for the ranges of 0.1 < U/k, J/k < 10. Pair- 
wise nonlinearity measured by U in the driven system is 
strong enough for the survival of pairwise entanglement 
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FIG. 2: (Color Online) Dependence of (a) ff (2) (0), (b) C (c) A12 for F/k = 0.1 and (d) 5 (2) (0), (e) ( (f) A12 for F/k = 1 with 
respect to dimensionless J/k and U / k in two-site KH system with two-photon exchange. £ is multiplied by 10 3 and A12 is 
multiplied by 10 5 for visibility. 



in the steady state. It is easier to violate spin squeezing 
inequality in the region of J/k < 1, where antibunching 
and hence particle like behavior is more significant. This 
is intuitively expected as wave like character enforced by 
the delocalizing photon hopping makes spin squeezing 
and particle entanglement more difficult, while nonlinear 
interaction favors pairwise particle correlations. 

Mode entanglement parameter A12 behavior with J 
and U, shown in Fig. 1(c) is opposite to that of spin 
squeezing. Lowest and highest values of A12 emerge in 
the antibunching and bunching regions of Fig. 1(a) re- 



spectively. Localization of photons in the cavities favor 
pairwise particle entanglement; while wave like behavior 
associated with coherent derealization by the photon ex- 
change favors mode correlations. Though A12 approaches 
toward positive values as J / k increase; mode correlations 
remained in the steady state are not strong enough for 
entanglement for the range of J/k we considered. In com- 
parison to Fig. 1(b) we see that particle correlations are 



more robust against cavity loss. 



Fig. 1(d) Fig. 1(e) and Fig. 1(f) depict quantum co- 
herence, spin squeezing and mode entanglement for the 
case of strong drive of F/2ir = 0.4 MHz, respectively. Co- 
herent effects of the drive is harmful for the particle like 
character and the range of g^ 2 '(0) is shifted away from 
and becomes 0.00486 < £ (2) (0) < 1.0321. Stronger non- 
linearity is needed for antibunching and thus the region 
of antibunching is narrowed in the U/k axis relative to 
that in Fig. 1(a) Spin squeezing and pairwise entangle- 
ment becomes harder to establish relative to weak drive 
case. The maximum value of £ reduces about an order 
of magnitude when the coherent drive is increased an or- 



der of magnitude. Spin squeezing is almost absent in 
the bunching region. On the other hand, mode entangle- 
ment is approximately established in the same bunching 
region. 

It should be stressed that the single photon exchange 
model gives clearly distinct roles to the interactions. 
Nonlinearity is localizing and a single-axis twisting spin 
squeezing type of interaction. Photon hopping is delo- 
calizing and a mode correlating interaction. Two photon 
hopping however has particle and delocalizing character. 
It is delocalizing due to photon exchange effect, but also 
establishes pairwise correlations by a two axis twisting 
interaction. We will point out emergence of distinct re- 
lation between second order coherence and entanglement 
in the following subsection. 



B. Two-Photon Exchange 

In Fig. [21 we show our results for two-photon exchange 
interaction. Fig. 2(a) Fig. 2(b) and Fig. |2(c) depict 
quantum coherence, spin squeezing and mode entangle- 
ment for F/2ir = 0.04 MHz, respectively. In contrast 
to single photon exchange, there is no coherent bunch- 
ing region in Fig. 2(a) and g( 2 '(0) varies in the range of 
0.00723 < g {2) {0) < 0.7944. The antibunching is strong 
in both delocalizing photon hopping and local cavity non- 
linearity dominated regimes. Relatively weak antibunch- 
ing occurs in J/k, U/k < 1 region where coherent drive 
is comparable to these local and nonlocal pairwise inter- 
actions. 

In the two photon exchange case, analogous to LMG 
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model, there are both single axis and two axis twisting 
routes to spin squeezing. When J/k > 1 the former 
becomes the major route leading to spin squeezing; while 
when U/k > 1 the latter becomes the main interaction 
causing spin squeezing. Accordingly strong violation of 
spin squeezing occurs over the entire domain of J/k, U /k. 
Similar level of violation can only be found in the narrow 
strip of weak J/k in Fig. 1(b) The uniformity of the 



pairwise entanglement witness parameter with U which 
is present in Fig. 1(b) is no longer present in Fig. |2(b)| 
Variation of spin squeezing and antibunching strength 
with hopping and nonlinearity becomes identical for the 
case of two photon exchange. 

Complementary character of mode and pairwise entan- 
glement is again revealed in comparison of Fig. 2(c) and 



Fig. 2(b) Regions of minimum violation of spin squeez- 
ing inequality coincides with the regions where mode en- 
tanglement parameter is maximum. Mode correlations 
are much weaker in the case of two photon exchange as 



can be seen by comparing the maximum of C in Fig. 2(c) 
to that of Fig. 1(c) Variation of g^{0) with J/k and 



U/ k in comparison to that of A12 turns out to be of recip- 
rocal nature for two photon exchange case. The largest 
mode entanglement parameter regions are the same with 
strongest antibunching regions. This reflects the particle 
like character of derealization induced by two photon 

hoppi ng. 

Fig. |2(d)| Fig. 2(e) and Fig. 2(f) show quantum coher- 
ence, spin squeezing and mode entanglement for F/2tt = 
0.4 MHz, respectively. Strong drive helps to bring coher- 
ence into the system and now the maximum g^ 2 '{0) ~ 1. 
The separation of two strong antibunching regions by a 
weak antibunching region around the line J — U becomes 
more visible. Similar lines divides the spin squeezing and 
mode entanglement figures into regions with same statis- 
tical character as well. Strong drive works against spin 
squeezing and pairwise entanglement; while it helps to 
establish mode correlations. The correlations are still 
not entangled in Fig. 2(f)| and pairwise entanglement, 
though becomes harder to establish, survives in Fig. 2(e) 



IV. CONCLUSION 

We considered single photon and two photon exchange 
coupled cavities with Kerr nonlinearity under drive and 
dissipation. We examined the steady state second order 
coherence, spin squeezing and pairwise entanglement as 
well as mode entanglement. 

Complementary character of pairwise and mode entan- 
glement is revealed. Drive strength can be used to tune 
between mode and particle correlations. Relations be- 
tween localization and bunching as well as derealization 
and antibunching is found to be translated to pairwise 
and mode entanglement in the case of single photon ex- 
change. Due to its nonlocal but pairwise interaction na- 
ture, such direct relations becomes more subtle in the 
case of two photon exchange, where derealization has 
particle character. Equivalent behavior with nonlinear- 
ity and photon hopping is found between coherence and 
spin squeezing. Pairwise entanglement and strong anti- 
bunching ranges are enhanced towards the strong photon 
exchange regime with the two photon coupling. In addi- 
tion to identification of coherent and correlated photons, 
second order coherence can discriminate the type of these 
correlations. 

Our results illuminates the subtle relations between lo- 
cality, coherence and quantum correlations that could be 
exploited for synthesis of many body quantum entangled 
states in coupled cavities with designed short and long 
range interactions. 
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